NASA TECHNICAL 






LOAN COPY: RETURN l£§S 
AFWL (WLIL-2) 
KIRTLAND AFB, N 


COUPLED VIBRATION AND 
DISSOCIATION RELAXATION 
BEHIND STRONG SHOCK WAVES 
IN CARBON DIOXIDE 

by Franz J. Hindelang 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. • FEBRUARY 1967 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 


0DbaM7fl 
NASA TR R-253 


COUPLED VIBRATION AND DISSOCIATION RELAXATION BEHIND 

STRONG SHOCK WAVES IN CARBON DIOXIDE 

By Franz J. Hindelang 

Ames Research Center 
Moffett Field, Calif. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $2.00 



COUPLED VIBRATION AND DISSOCIATION RELAXATION BEHIND 


STRONG SHOCK WAVES IN CARBON DIOXIDE 

By Franz J. Hindelang 
Ames Research. Center 


SUMMARY 


The harmonic oscillator rigid-rotator model has been used to calculate 
the relaxation region behind a shock wave in carbon dioxide. Finite relaxa- 
tion rates for the three different vibrational modes and two dissociation 
reactions are included. Models for the coupling between the vibrational relax- 
ation and the dissociation process are based on the assumption that dissocia- 
tion can proceed from any vibrational level with equal probability. Two 
different models for the vibrational excitation have been examined. Solutions 
have been obtained for the interdependent fluid-flow, chemical rate, and 
vibrational relaxation-rate equations incorporating estimated rate coeffi- 
cients. Results are presented in the form of flow-field profiles for density, 
pressure, translational and vibrational temperatures, and species concentra- 
tions. The effects of vibrational excitation, vibration-dissociation coupling, 
and energy exchange between the vibrational modes are investigated. The 
effect of vibrational relaxation and vibration-dissociation coupling is much 
stronger in CO 2 with three different vibrational modes than in diatomic gases 
with only a single mode. The results of this study show that the effect of 
coupled vibrational relaxation and dissociation can sometimes alter the flow- 
field profiles by a factor of 2 compared to similar calculations without such 
coupling. For vibrational relaxation the results indicate that shock-wave 
profiles depend primarily on the rate at which the translational energy is fed 
into internal modes and not so strongly on the energy distribution among the 
modes . 


INTRODUCTION 


The atmospheres of Mars and Venus are composed of a mixture of gases of 
which C0 2 is a known constituent. The effect of even a small amount of C0 2 on 
aerodynamic heating and flow about a hypervelocity entry body can be signifi- 
cant because polyatomic molecules like CO 2 can be excited to much higher 
vibrational energies than diatomic molecules at the same temperature and will 
dissociate at lower temperatures. To predict the properties of the flow field 
about a vehicle entering an atmosphere containing C0 2 and the aerodynamic 
heating it will experience, it is necessary to understand the high temperature 
chemical kinetics of C0 2 , particularly, the dissociation and vibrational 
relaxation processes and the manner in which they couple. 

Estimates of the flow field behind strong shock waves in CO^ with 
dissociation relaxation have been obtained by Howe, Viegas, and Sheaffer in 
reference 1, and with an interdependent dissociation and ionization relaxation 



process by Howe and Sheaffer in reference 2 . These references used a 
collision theory to predict the reaction-rate coefficients. This paper is 
basically an extension of the work of Howe et al. to include the effects of 
vibrational relaxation and vibration-dissociation coupling in C 0 2 . These addi- 
tional effects are important for calculating more accurately the flow field 
quantities and their gradients immediately behind the shock front, where C0 2 
dissociation takes place. For example, the translational temperature immedi- 
ately behind the shock front is about twice as high when the vibrational 
energy is unexcited compared to the same shock wave with completely 
equilibrated vibrational energy. 

Particular emphasis is given in this study to the nonequilibrium of three 
different vibrational modes. This complicates the analysis significantly in 
comparison to shock wave calculations for diatomic gases (refs. 3 “ 5 )- On the 
other hand, the relatively simple, harmonic oscillator rigid-rotator model is 
used to describe the coupling between vibration and dissociation, thus avoid- 
ing more complicated models. It is assumed that dissociation can proceed from 
any vibrational level with equal probability if a collision provides suffici- 
ent energy. The population of the vibrational levels is represented by 
Boltzmann distributions with effective vibrational temperatures. Numerical 
values for the dissociation rates of C 0 2 and CO are taken from reference 1 . 
Vibrational relaxation times at high temperatures are extrapolated from the 
available experimental results. 

The research reported in this paper was conducted while the author was 
pursuing a NAS-NRC Postdoctoral Resident Research Associateship supported by 
the National Aeronautics and Space Administration. The support from the 
National Academy of Sciences and NASA is gratefully acknowledged. 


SYMBOLS 


Units are generally on a per mass basis (superscript A indicates per 
mole basis) . 

B c coefficient in backward reaction-rate expression (eq. (27)) for 
reaction r 

Bp coefficient in forward reaction-rate expression (eq. (26)) for 
reaction r 

Cp specific heat at constant pressure 

Cy specific heat at constant volume 

Dj_ dissociation energy of species i 

E internal energy of gas mixture 

E Cr defined by equation (27) 
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Eg activation energy in forward reaction-rate coefficient (eq_. ( 26 )) for 

r reaction r 

internal energy of species i 

e vibrational energy of gas mixture 

e± m _brational energy of species i in mode m 

^ijp average vibrational energy in mode m of species i (lost by 
1 dissociation) 

f number of active degrees of freedom in C0 2 

h static enthalpy or Planck 1 s constant 

Kc r equilibrium coefficient for reaction r 

K m defined by equation (Al) 

k Boltzmann* s constant 

kb r backward reaction-rate coefficient for reaction r 

k-p forward reaction-rate coefficient for reaction r 
± r 

M Mach number or typical collision partner 

N number of bimolecular collisions per second to which a C0 2 molecule is 

subjected 

nj_ moles of species i per unit mass of mixture 
P defined by equation (Bl) 

P m transition probability for mode m 

p static pressure 

R universal gas constant (per mole basis) 

T absolute temperature 

273° K 

T' defined by equation (B4) 

u velocity in x direction 

V m total vibrational levels in mode m of C0 2 

v m vibrational level in mode m of C0 2 
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a chemical species 
distance behind shock front 
exponent in equation (27) 

forward stoichiometric coefficient and exponent in equation (9) and (10) 

backward stoichiometric coefficient 

Cp/c v , ratio of specific heats 

Pg/p^, density ratio across shock front 

defined by equation (31) 

hv/k, characteristic vibrational temperature 
vibrational frequency of one mode of C0 2 
mass density of gas mixture 
standard or sea-level density, 1.225 kg/m 3 
relaxation time for mode m of C0 2 

relaxation time for energy exchange between the valence and bending 
modes of C0 2 


Subscripts 

collision 

dissociation 

equilibrium 

species i, i = 1, 2, 3, k corresponds to C0 2 , CO, C, and 0, respectively 

vibrational mode m or m r (l, valence; 2 and 3, bending; k y asymmetric 
stretching mode) 

rth reaction as shown in equations (2^) and (25) 

conditions immediately behind shock front (only translation and 
rotation in equilibrium) 

vibrational level in mode m 

conditions ahead of shock 


k 



Superscripts 


per mol basis 
mean value 


ANALYSIS AND SOLUTION 


Before beginning the detailed description of the analysis of this paper, 
a qualitative discussion of the problem and the assumptions involved in the 
analysis are given. When finite exchange rates for the internal energies of 
the molecules are considered, the shock can no longer be assumed to be a 
vanishingly thin discontinuity. In fact, its processes can extend consider- 
ably in the direction of the flow. A schematic temperature profile of the 
flow field to be analyzed in the following sections is shown in sketch (a), as 
seen by an observer stationary with the shock front. 



In a shock wave, the translational and rotational energy modes become 
excited first. Only a few collisions fully excite these modes, so they are 
usually taken to be fully excited within a narrow region defined as the shock 
front. Vibrational excitation requires many more collisions and, at suffici- 
ently low pressures, a finite relaxation zone is evident. Dissociation is 
coupled with vibrational excitation since it can occur from all vibrational 
levels. At low shock speeds and their associated low translational tempera- 
tures, dissociation does not proceed appreciably until vibrational excitation 
is near equilibrium, since here dissociation proceeds predominantly from 
higher vibrational levels. At high shock speeds, that is, high translational 
temperatures, more energetic collisions occur causing more dissociation from 
lower vibrational levels. Consequently, vibrational relaxation and dissocia- 
tion proceed at comparable rates. Ionization could, of course, require an 


5 



additional relaxation length, hut it is not considered in this paper, since 
at the highest shock speed for which flow field profiles have been calculated, 
u^ = 9 km/s, the equilibrium electron concentration is only about 2 percent. 

Although there are many possible reactions for the dissociation of carbon 
dioxide, in order to simplify the thermodynamic and kinetic models, only the 
following reactions are considered in this study: 

C0 2 2 GO + 0 
CO ^ C + 0 

Vibrational relaxation and vibration-dissociation coupling are considered for 
the C0 2 molecule only, and the relaxation region is analyzed to the point 
where alnost all the C0 2 is dissociated. 

The coupling between vibrational and dissociation-relaxation processes in 
diatomic gases has been investigated by various authors (refs. 3-5). In this 
report, flow fields are investigated with vibration-dissociation coupling in 
the polyatomic molecule C0 2 . The model used is similar to that developed for 
diatomic gases, but several further assumptions are made: 

(l) A harmonic oscillator model with a constant moment of inertia for 
rotation (rigid rotor) is assumed for the vibrational energy of C0 2 . A linear 
molecule like C0 2 has four different normal modes, two of which are degenerate, 
thus giving three different vibrational frequencies, as shown in sketch (b) . 
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Sketch (b) 


(2) The transition between levels in each mode as well as energy exchange 
between modes is described by an appropriate relaxation equation with a char- 
acteristic relaxation time instead of many transition probabilities. A com- 
plete mathematical description of the excitation of the vibrational states and 
the dissociation from these states requires knowledge of all transition proba- 
bilities between all eigenstates, including the continuum dissociation state. 
Since there are no reliable data for the transition probabilities from theory 
or experiment, such an analysis would involve too many new unknowns, and the 
above-mentioned simplification was made to obtain a tractable form for 
solution. 

(3) Only the ground electronic state was assumed to occur for C0 2 . 

(4) The reaction products C and 0 were taken to be calorically perfect 
(c v = const), which simplifies the problem. This is justified because the 
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prime objective of the study is to show the effects of vibration-dissociation 
coupling in C0 2 on the flow field where C0 2 is present with little concern 
about processes far downstream. 

( 5 ) The vibrational energy of CO was assumed to be in equilibrium with 
the translational temperature. This raises the important question of the 
vibrational energy of a diatomic molecule formed from the dissociation of a 
polyatomic molecule. It can be assumed that the CO molecule would carry away 
at least part of the vibrational energy of the C0 2 molecules. The initial 
state of the CO molecule, as well as subsequent vibration-dissociation 
coupling in CO, requires further study. 

The mathematical problem is to solve the conservation equations and the 
equations of state with the appropriate rate equations. This is outlined in 
the next section. 


Basic Equations 

The flow field behind the shock wave can be described by the equations of 
conservation of mass, momentum, and energy, the appropriate equations of state, 
and the chemical and vibrational-relaxation equations. The solution of this 
system of equations determines the thermodynamic variables, the velocity of 
the fluid flow, the species concentrations, and the vibrational energies of 
the C0 2 molecules. 

When transport phenomena are neglected, the conservation equations are: 

P u = Poo U oo = const (l) 

( 2 ) 

( 3 ) 

The model of an ideal gas gives a good approximation for the thermal 
equation of state of gases at high temperatures: 


du dp 

pu — — + — — = 0 
dx dx 


du , dh ^ 
dx dx 


P = pRT 2^ n ± 
i 

The internal energy of a mixture of ideal gases is: 

E 


(b) 


" I ® irii 


( 5 ) 


where 


7 



“ c vi^ + 


( 6 ) 



m 



The term c v is the specific heat for changes in translational and rotational 
energy only, and is assumed constant at its fully excited classical value (3/2 
R for atoms, 5/2 R for linear molecules). The term is the contribu- 

tion to the vibrational energy of the species i from the vibrational mode m. 
It is obtained by summing over the population of all levels in mode m (see 
appendix B) . 

Since electronic excitation is ignored, the enthalpy depends only on 
pressure, density, chemical composition, and the vibrational state of the 
molecules. It is given by 

h h-jHi (e im + Sp.TiK (7) 

i i,m 

where 

g pi = R + c Vi 


To complete the set of equations describing the flow field behind the shock, 
the time dependent vibrational and chemical rate equations are included. The 
vibrational rate equation has the form: 


d§i 

1 m 


dx 





(8) 


The first term on the right-hand side of equation (8) is the usual vibra- 
tional relaxation term for collision-induced energy exchange between transla- 
tional and vibrational degrees of freedom and between different vibrational 
modes. The second and third terms account for vibration and dissociation 
coupling. The term Sj_ m - §i m is the net average vibrational energy of 
species i in mode m lost by dissociation. It is obtained by summing over 
all levels, v m (see eq. ( 29 ))* The term e± m accounts for the change in the 
number of oscillators during dissociation. It is assumed that a molecule can 
dissociate with equal probability from any vibrational level, if a collision 
provides enough energy. Since higher vibrational levels require less energy 
exchange during a collision, dissociation proceeds predominantly from higher 
vibrational levels. This means that a steady loss in vibrational energy 
occurs as long as dissociation proceeds, since > §j_ . 

Similarly, one can expect recombination to increase vibrational energy. 

In this study, the effect of recombination on the vibrational energy is 
neglected since for flow fields where dissociation is the predominant reaction, 
such as the dissociation of C0 2 behind strong shock waves, recombination con- 
tributes very little to the vibrational energy. Indeed, Treanor and Marrone 
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(ref. k) found that for shock waves in diatomic gases, like N 2 and 0 ^, the 
contribution of recombination is negligibly small until conditions extremely 
close to full equilibrium are reached. Therefore, we neglect the effect of 
recombination on vibration. 


For a collision controlled reaction r such as 

kf T 


y a r j i^i = V Pr.iXi 


(9) 


r u r 1 

The chemical reaction-rate equation can be written in the usual form: 


r 


(io> 


The dissociation rate coefficient kf is a function of the vibrational 
state of the dissociating molecules and of the translational temperature. 

Thus, the interaction between vibration and dissociation enters the system of 
equations through kf as well as the vibrational relaxation equations. The 
coefficients k^ are obtained by relating the forward and backward rates to 
the equilibrium constant. This is permissible for a simple harmonic oscil- 
lator even for the conditions of -vibrational nonequilibrium (ref. 6 ). For 
polyatomic molecules, it is the only available approximation to obtain recom- 
bination rates for vibrational nonequilibrium. It should be pointed out that 
recombination rates are unimportant as long as the equilibrium concentration 
of the recombining species is negligible throughout the flow field as is true 
in the study for all but the low shock-speed case (u^ = 2.5 km/s). 

The boundary conditions of the flow and chemical-rate equations at the 
beginning of the relaxation zone (x = 0 ) are obtained from the normal shock 
equations of an ideal gas by assuming that within the shock front no chemical 
reactions occur and the vibrations are unexcited. The ratio of specific heats 
7 is therefore constant (7 = 7/5 for linear molecules). The boundary 
conditions have the form: 


with 


P = P s = P„ 


2^ - (y - l) 

7 + 1 


= UaVPoo/yPco 

. = . = _ (z + DHx> 2 

P P s p w (7 + l)!^ 2 + 2 


(11) 


( 12 ) 


h — h s 


, 1 
— h + “ 

00 2 


(Uco 2 “ U s 2 ) 


(13) 
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( 14 ) 


^im ^ijii^oo 


ni = n lc 


n.2 = n 3 = n.4 = 0 


(15) 

(16) 


where ni, n 2 , n 3 , n4 refer to the concentrations of C0 2 , CO, C, 0, 
respectively. 

The equations and appropriate boundary conditions for describing the 
flow field behind the shock front is thus complete. 

In the following sections, details of the vibrational relaxation, the 
chemical reactions, and the interaction between vibration and dissociation are 
discussed. 


Vibrational Relaxation 


To evaluate the terms (dei m /dx) c in equation (8), a harmonic oscillator 
model is used to describe the vibrations. The vibrational relaxation is 
limited to the C0 2 molecules. 

Consistent with the model of a harmonic oscillator, the CO2 molecule is 
taken to have four vibrational degrees of freedom and, accordingly, four nor- 
mal modes of vibration, as shown in sketch (b). (A linear molecule with N 
atoms has three N-5 degrees of freedom.) Two of these degrees of freedom are 
degenerate, however, leaving three different fundamental frequencies of vibra- 
tion. The vibrational frequencies in this study are taken from reference 7* 
For CO2 molecules, an accidental degeneracy occurs between the valence (sym- 
metric stretching mode) and bending mode (vi « 2V2) • An energy exchange 
between these modes is termed resonance exchange, since no exchange with 
translational energy is involved. If the terms are retained for direct exci- 
tation of the symmetric and asymmetric stretching modes, which were omitted in 
reference 8, then the equations for vibrational relaxation can be written as: 

0 Er) * u [ ^ i(T) - + t«ii(T a ) - indnij (17) 


where §n(T) is the vibrational energy of C0 2 , mode 1 in equilibrium with 
temperature T . 


(%r) c = ^(T) - §12(12)] - ^ [Sii(T 2 ) - §n(Ti)]j 

f&ia\ = /de 12 \ 

V J c ~ \ <3x ) c 


(18) 

(19) 


10 



( 20 ) 


0Gr) c - k [ ^ (T) - ^< T ‘> '} 


Symmetric stretching mode 
(valence mode) 

m = 1 

0i = 1920 ° K 

Bending mode 

fm = 2 
W = 3 

02 = 03 = 960 0 K 

Asymmetric stretching 
mode 

m = 4 

04 = 3380° K 

The relaxation times. 

r , can 

be expressed in terms of transition 


probabilities (ref. 8) 


i = NP^f°(l - e- 0 «i/ T m) 
T m 


( 21 ) 


This equation is used to compute the flow-field profiles. The term P m is 
the probability that during a collision the vibrational mode m is de -excited 
from the first to the zeroth quantum state. Experimental results at low 
temperatures (T < Q m , see appendix A) are usually referred to an approximate 
form for T m . 



i-*o 

m 


( 22 ) 


The relaxation time for the energy exchange between the modes is: 


^ = 4P^Ng 12 (T 2 )[£ 12 (T 2 ) - en(T 2 ) ] -1 ( 23 ) 

T 12 

o-^i 

The term P g _ >Q is the probability that during a collision the valence 
mode is excited from the zeroth to the first quantum state while the bending 
mode is de -excited from the second to the zeroth state and is proportional 
to the translational temperature T (see ref. 8) . 

Values for the various vibrational relaxation times used in this paper 
are given in appendix A , where the experimental work on vibrational relaxation 
times is reviewed. It is believed that the lack of accurate relaxation times 
can be tolerated in this paper , since the quantitative effects of coupling 
between vibration and dissociation can be found by a parametric study. 


Chemical Reactions 

The dissociation and recombination of CO 2 and its constituents are 
described by the following reactions: 
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1: 

C0 2 + M ^ CO + 0 + M 

(24) 

2: 

CO + M ? C + 0 + M 

(25) 


where M is an arbitrary collision partner. The formation of 0 2 due to 
recombination has been omitted, since its contribution to the chemical behav- 
ior of C0 2 is negligible (see ref. 2). The rate coefficients of the forward 
reactions can be expressed in the Arrhenius form: 


a fr -E fr /RT 
kf r = B f r T e 


( 26 ) 


The values for the coefficients Bf r and are taken from reference 1 

and are listed in appendix A. 

The rates of the backward reactions, which are important only near 
thermodynamic equilibrium, can be obtained from the following expression, with 
the notation of reference 1: 


kf r / ac. 

kv. = p ? — = kf B c T 
D r K n u r 


A iA 

r - E c r /R' 
e 


-l 


(27) 


where K Cr , the equilibrium constant for the reaction, is a function of the 

1 /rT 

translational temperature. In the next section e r and e" m will be 

calculated as functions of the vibrational state of the molecules. 


Evaluation of e and e m 

The interaction between vibrational relaxation and dissociation enters 
the system of equations in two places, one in the forward dissociation rate 
expression (eq. ( 26 )), where the effect of vibrations upon dissociation enters 
and the other in the vibrational relaxation expression (eq. (8)), where the 
effect of dissociation upon vibration enters by draining vibrational energy. 
The effects of the interaction are concentrated in two terms, e“^f/^ and e m . 
The model assumed that a molecule can dissociate from any vibrational energy 
level with equal probability, if during a collision the relative line-of- 
center kinetic energy exceeds D x - e. Out of all molecules, a fraction pro- 
portional to P will have the vibrational energy e. The probability 
that a molecule with vibrational energy e dissociates during a collision 
is proportional to the probability that the line -of -center kinetic energy is 
equal or greater than D x - e . For a Boltzmann distribution in the transla- 
tional degrees of freedom, this probability is e”^i“ e VRT. Without vibration- 
dissociation coupling, the energy exchange during a collision must exceed D x , 
yielding the probability factor e -Di/BT # This factor would appear unchanged 
in the expression for the forward reaction rate (without coupling) since it is 
the same for all molecules. With vibration-dissociation coupling one has to 
account for the different vibrational energies of the molecules by summing 
over all energy levels: 
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( 28 ) 


& 


^ ^ p e "( I) l/ R-V l 0 l - ' V '202-V303-V404)/T 

( -if /ftp vi=o v 4 — o 




Vi=0 v 4 =o 


The energy removed from mode m with each dissociation is v m 0 m . If dissoci- 
ation occurs, then the probability^that it resulted from the levels Vi, v^^ 

v 4 is proportional to Pe Vl ^ 1 V 40 4 )/T^ ip^e average energy lost 

by dissociation from mode m is obtained again by summing over all energy 
levels. Equations (28) and ( 29 ) are thus obtained by applying the method of 
Treanor and Marrone (ref. 4) to polyatomic molecules. 


e 


m “ 



V 1-0 v 4 =o 



Vi=0 v 4 =o 


P - ( D i/R Vi0 1 - V202 “V303 -■ v 4 8 4 )/T 

e v mym 


p e -(fil/P-Vi0i-V202“‘V303-V404)/T 


( 29 ) 


with 


p _ e "[ (vi 0 i/Ti)+. . .+( v404/T 4 ) ] 

To evaluate equations ( 28 ) and (29),- it is necessary to count the 
vibrational levels v. This imposes a restriction on the total vibrational 
energy of a molecule. The total vibrational energy of a molecule cannot 
exceed the dissociation energy. This is equivalent to the cutoff harmonic 
oscillator model for diatomic molecules. The classical vibrational energy of 
C0 2 is a linear function of the vibrational temperature without an upper limit; 
e/R = 4T V . At temperatures above 16,000° K the classical model has no physi- 
cal meaning because it predicts vibrational energies in excess of the dissoci- 
ation energy. Two different models that limit the amount of vibrational 
energy to the dissociation energy are investigated in this study. The first 
is based on the postulate that the vibrational energy of each mode cannot 
exceed Di/4 and on the fact that spectroscopic data (ref. 7 ) do not show high 
vibrational energy levels. We conclude from this that it is unlikely that 
all vibrational energy will be found in one mode only. The second model Is 
based on the postulate that the total vibrational energy cannot exceed the 
dissociation energy Di but without restricting the distribution of the 
vibrational energy among the different modes. This means that all combina- 
tions of ei ... e 4 are allowed as long as e x + . . . + e 4 < D x » For a 
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3 


4 


v, = l 



molecule with two vibrational inodes, 
the possible vibrational energies are 
shown in sketch (c) with the same 
vibrational frequencies for both modes. 

It is instructive to discuss the 
limiting value of the vibrational 
energy for a mole of gas and to compare 
this with the classical value. Since 
the energy of each molecule is limited, 
the vibrational energy for a mole of 
C 0 2 at high temperature cannot exceed 
upper limits of approximately l/2 and 
4-/5 respectively, for the two models. 
These numbers are obtained by replacing 
the summations over all vibrational 
levels by an integration (see 
appendix B) . 


It is believed that the vibrational cutoff of a C 0 2 molecule will be 
somewhere between the two cases investigated. The mathematical formulation 
of both models is given in appendix B. 


Method of Solution 


Two different solutions for the nonequilibrium flow field have been 
obtained and are referred to as the complete and the approximate solution. 


The approximate solution is based on a partial uncoupling of the flow and 
the rate equations by an extension of a method given by Bethe and Teller 
(ref. 9 ). This technique proceeds by solving the rate equation (eq. (lO)) to 
find the time element during which a given fraction of the total C 0 2 disso- 
ciates with the assumption that the flow-field quantities T, p, and p are 
constant. The change in vibrational energy during this time element is given 
by the relaxation equations (17) through ( 20 ). From the chemical composition 
and vibrational energy at a particular time in the flow field, the time inde- 
pendent equations (l) through (7), are used in the following form to solve for 
the flow-field quantities. 


u 

a 


gT] 


/(gTl) 2 - (2t) - 1) 



2t] - 1 


(30) 


with 



i 


(31) 


14 



( 32 ) 


l + yMoo g 
yM 

' oa 


p 



Poo 

U 


JL - 

Poo 

7 M oo 1 

f u\ 

" 0 


(33) 

(3*0 


The next integration step starts with the new flow-field quantities. Results 
of the approximate solution were used to estimate the step size required for 
the Runge-Kutta integration process incorporated in the complete solution. 


The complete solution is essentially a numerical integration of the flow- 
field equations, the chemical reaction rate, and the vibrational relaxation 
equations subject to boundary conditions (ll) to (l 6 ) . The basic differential 
flow relation for numerical integration is obtained by combining equations (l), 
( 2 ), ( 3 ), (k), and ( 7 ), making x the independent variable (ref. 10 ) . The 
result is given by equation ( 35 )* 


du 

dx 




i 


u 



n i c Pi 



dn-^ 


~Pi dx 



i 


n i c Pi 





n 


dei 

idx 7 



(35) 


In equation (35) is the specific heat for changes of translational and 

rotational energy only and is assumed to be constant at its fully excited 
classical value (5/2 R for atoms, 7/2 R for diatomic and linear triatomic 
molecules) . 


Both methods of solution gave, within numerical errors, the same results. 


DISCUSSION OF RESULTS 

Solutions of the nonequilibrium flow field for shock waves in pure C0 2 
were obtained and are discussed in this section. Before the results are 
examined, however, a short discussion of vibrational relaxation times and a 
dissociation relaxation time is presented to provide the reader with a better 
understanding of the shock speeds at which vibration dissociation coupling is 
important. The dissociation relaxation time is defined by an equation similar 
to the simple vibrational relaxation equation, that is. 
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( 36 ) 


dni (n x ) e - n x 

Immediately behind the shock front, recombination can he neglected, and 
(n x ) e = 0, giving 

RT 

T a? = % 

A graphical comparison of the relaxation times is shown in figure 1 as a 
function of the translational temperature immediately behind the shock wave. 

The experimental data at low shock 
speeds were extrapolated to high tem- 
peratures. It is estimated that these 
extrapolated vibrational relaxation 
times are too high because transitions 
between more than one vibrational level 
can be induced by collisions with high 
line -of -center relative energy. Fur- 
thermore, near collisions of molecules 
may also induce vibrational transitions. 
Both effects are not included in the 
model used to describe the vibrational 
relaxation. The extrapolation cannot 
give reliable relaxation times at high 
temperatures. This study is therefore 
parametric. The extrapolated relaxa- 
tion times are used as base values, and 
have been multiplied by constant fac- 
tors so as to show the effects of vibra- 
tional equilibrium, vibrational 
relaxation, and unexcited vibrations on 
the flow field at various shock speeds. 

Results have been obtained for three different shock speeds. For low 
shock speeds, u^ < 3 km/s (T s < 7000° K) , the dissociation is slower than the 
vibrational excitation rate, thereby allowing vibrational equilibrium to be 
achieved before dissociation becomes appreciable. At intermediate shock 
speeds between 3 and 6 km/s, the coupling between vibration and dissociation 
is expected to be more important. At high shock speeds, u^o > 6 km/s, addi- 
tional processes like CO dissociation and ionization become more and more 
important. In this case, the region of C0 2 dissociation is only a small part 
of the whole relaxation region. 

Typical examples of flow-field solutions are discussed in the following 
sections . 
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Low Shock Speed 

The flow-field profiles for low shock speeds, u s < 3 km/s, are character- 
ized by insignificant amounts of dissociation. At these low shock speeds. 


vibrational relaxation is predominant. 



O 05 I 15 2 

Distance behind shock, x, cm 


(a) Density ratio and pressure ratio. 



o 05 1 15 2 


□ ■stance behind shock, x, cm 

(b) Translational and vibrational 
temperature . 

Figure 2.- Nonequilibrium flow-field 

profiles. Effect of interaction between 
the valence and the bending mode for low 
shock speeds (u = 2.5, P^/p^ = 10 “ 2 , 

Too = 273° K) . 


A typical solution for the flow-field 
profiles with vibrational relaxation is 
shown in figure 2 for the shock condi- 
tions Ucq = 2.5 km/s, density 
Poo /p o = 10-2, and ^ = 273 ° K. The 

results are given for varying inter- 
action between the valence and bending 
modes as follows: (a) fast energy 

exchange between the two modes, 0, 

(b) energy exchange between the two 
modes with a relaxation time ti 2 ~t 2 , 
and (c) no energy exchange between the 
two modes, t X2. °°* The relaxation times 

for the valence, bending, and asymmetric 
stretching modes, ti, t 2 , t 4 , respec- 
tively, are given in appendix A. In this 
model t 2 is about 1 order of magni- 
tude less than ti or t 4 which means 
that the bending mode is excited first. 

Pressure and density profiles are 
shown in figure 2(a). The pressure is 
relatively insensitive to the overall 
relaxation process and even more insen- 
sitive to the energy exchange between 
the valence and bending modes. The 
density is quite sensitive to the over- 
all relaxation process and only 
slightly sensitive to energy exchange. 

The translational and various 
vibrational temperatures are shown in 
figure 2(b). The temperature of the 
valence and bending mode is quite sen- 
sitive to the energy exchange between 
the modes. For no energy exchange, the 
temperature of the bending mode over- 
shoots the translational temperature. 
The translational temperature is also 
slightly lower when an equilibrium 
energy exchange between the valence and 
bending mode occurs. The variation in 
translational temperature is the main 
driving force for the variation in 
density shown earlier, since the pres- 
sure is virtually unaffected. 
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The variation of the vibrational 
energies in the valence, m =• 1, bending, 
m = 2, 3, and asymmetric stretching 
mode, m = 4, is shown in figure 2(c) to 
follow the trends indicated *by the 
vibrational temperatures. It should be 
noted that the bending mode has the 
highest energy and the valence and 
bending modes relax faster than the 
asymmetric stretching mode. All the 
curves were terminated at x = 0.2 cm 
for convenience and should not be taken 
to imply that vibrational equilibrium 
occurs in this distance. The equilib- 
rium values are noted on the figures. 


(c) Vibrational energies. 

Figure 2.- Concluded. 

Intermediate Shock Speed 

Profiles for an intermediate shock speed, u s = 5 km/s, p 00 /p Q = 10 4 , and 
T^ = 273 ° K, are shown in figure 3. These curves were calculated for vibration- 
dissociation coupling with the vibrations (a) frozen at the free-stream vibra- 
tional energy with dissociation proceeding predominately from the ground state 
of vibration, (b) relaxing at a finite rate (the relaxation times were multi- 
plied by 10 " 1 to show the coupling), or (c) relaxing at an infinite rate to 
keep the vibrational energy in equilibrium with the translational temperature. 
The effect of different limiting energies is also shown; the limiting cases 
being that the vibrational energy of a molecule cannot exceed D/k per mode 
(model I), and that the vibrational energy of the molecule cannot exceed D 
without further restriction (model II). 

The total vibrational energy as well as the vibrational energy per mode 
is shown in figure 3(a) for case (b). Again it is noted that the bending 
modes contain by far the largest portion of the total vibrational energy, 
especially when their degeneracy is considered. 

The translational temperature behind the shock wave is shown in 
figure 3(b). As expected, the combination of variables that gives the highest 
vibrational energy (t 0 and model II, i.e., vibrational energy of a molecule 
not exceeding D) also gives the lowest translational temperature. 

The translational temperature immediately behind the shock front with 
vibrations in equilibrium is lower than for frozen vibrations, see appendix C. 

The pressure is again shown in figure 3(b) to remain nearly constant 
behind the shock wave with the combined variables giving the highest vibra- 
tional energy as well as the highest pressures. Consequently, the density is 
also higher for higher vibrational energy. The density ratios for vibrations 
frozen or in equilibrium differ by as much as 100 percent (see fig. 3(c)). 
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The mole fractions of CO and C0 2 behind the shock wave are shown in 
figure 3(d). It is noted that dissociation of C0 2 proceeds much more rapidly 
when the vibrational energy is higher, since more molecules are available in 
the higher vibrational states and consequently dissociate more readily upon a 
collision with another particle. The numerical results for the two models do 
not differ too much even though the assumptions on vibrational cutoff are 
quite different. The assumptions affecting the rates of vibrational relaxa- 
tion (i.e., vibrations frozen, relaxing, or in equilibrium) have a greater 
effect on the flow field than the assumptions on vibrational energy limit. 
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(a) Vibrational energies (in the case 
of vibrational relaxation) . 







(c) Density ratio. (d) Species concentration. 

Figure 3 .- Effect of vibrational energy limit on nonequilibrium flow-field profiles for 
intermediate shock speeds (u^ = 5 km/s, P^/pQ = 10~S = 273° K) . I. Vibrational 

energy of a molecule not exceeding D/4 per mode (d/ 2 for the bending mode). 

II. Vibrational energy of a molecule not exceeding D. 
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High Shock Speed 


Flow-field profiles for a high shock speed, u^ = 9 km/s, P^/Pq = 10 “ 2 , 
and = 273° K, are presented in figure 4 for the limiting cases of frozen 
and equilibrium vibrations as well as for finite relaxation times (relaxation 
times multiplied by 10“ 3 ). Two different sets of data are given: one set 

with coupling between vibration and dissociation and with the vibrational 
energy per molecule limited to 1/4 D per mode, the other with no coupling and 
with the vibrational energy unrestricted, that is, the vibrations can be 


excited to the full classical value. 



(a) Temperature ratio. 



(c) Density ratio. 

Figure 4.- Nonequilibr lum flow-field 
profiles (u^ = 9 km/s, p OT /p = 10 
= 273 0 K) . Effect of interaction 
between vibration and dissociation. 


In this case it is assumed that 



(b) Pressure ratio. 


dissociation proceeds from the ground 
state of vibration. The results for 
vibrational equilibrium in this case 
are then identical with the earlier 
work of reference 1. 

The translational temperature 
behind the shock wave is shown in fig- 
ure 4(a). The temperature immediately 
behind the shock wave for vibrational 
energy frozen or relaxing is higher 
than when vibrations are in equilibrium 
and is, in fact, the ideal gas value 
with only translational and rotational 
equilibrium excitation. The difference 
between the temperature immediately 
behind the shock for vibrations in 
equilibrium with and without interac- 
tion shows the effect of limiting 
vibrational energy per molecule; the 
unlimited vibrational energy solution 
gives a lower translational temperature, 
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Distance behind shock, x, cm 

(d) Species concentration, interaction 
between vibrational and dissociation. 



Distance behind shock, x, cm 

(e) Species concentration, no interaction 
between vibration and dissociation. 


Figure 4.- Concluded. 

since more energy can "be fed into vibrations. The time required to obtain 
vibrational equilibrium is also shorter for the no interaction solutions. The 
behavior of the pressure and density behind the shock wave is shown in fig- 
ures 4(b) and 4(c) . A somewhat unexpected result for no interactions is that 
both pressure and density can have negative slopes. The solutions for 
vibration-dissociation coupling, however, show only positive pressure and den- 
sity profiles in the flow field. The mole fractions of C0 2 , CO, and 0 are 
shown in figures 4(d) and 4(e) . A significant effect of coupling is to 



Figure 5*- Nonequilibrium flow-field 
profiles (u w - 9 km/ s J Poo/p q = 1° 

T = 273° K) . Effect of large changes 
of the vibrational relaxation times. 
Temperature ratio and pressure ratio. 


increase the forward reaction rate, as 
can be seen by comparing the concentra- 
tion profiles with coupling, fig- 
ure 4(d), and without coupling, 
figure 4(e) . For the solutions without 
coupling, the concentrations are not 
very sensitive to vibrational relaxa- 
tion since vibrational excitation 
changes only the translational tempera- 
ture and not the activation energy Ep . 
Since the vibrational relaxation times 
used are extrapolated from low tempera- 
ture data, it is instructive to examine 
the sensitivity of the solutions to the 
magnitude of the relaxation times. A 
comparison of the temperature and pres- 
sure profiles behind the shock wave, 
with vibration-dissociation coupling, 
is shown in figure 5 for several values 
of vibrational relaxation times. It is 
seen that the vibrational relaxation 
time would have to be at least two 
orders of magnitude smaller than those 
assumed for this study before a signif- 
icant effect upon the flow field could 
be seen. 
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CONCLUDING REMARKS 


The coupling between vibration and dissociation relaxation behind shock 
waves in C0 2 has been studied for shock speeds from 2.5 to 9 km/s. Solutions 
of the fluid flow, chemical rate, and vibrational relaxation equations have 
been obtained, assuming frozen vibrations, vibrational relaxation, and vibra- 
tions in equilibrium both with and without coupling to the dissociation 
process. Results are presented as flow-field profiles for the flow variables 
and species concentrations. 

At shock speeds below 3 km/s, vibrational relaxation is the predominant 
relaxation process in the flow field, since dissociation is negligible. At 
higher shock speeds, the beginning of the relaxation region is governed by 
the interdependent vibrational relaxation and dissociation processes. The 
temperature immediately behind the shock front depends strongly on vibrational 
excitation. For frozen vibrations, the translational temperature at the shock 
front is more than twice as high as for instantaneous vibrational excitation 
to the full classical value. The parametric study shows that vibration- 
dissociation coupling in C 0 2 can alter the forward reaction rate considerably. 
Therefore, vibration-dissociation coupling must be considered when experi- 
mental rates are extrapolated over a wide range of temperatures. The esti- 
mates of rates in this study indicate that coupling becomes effective for 
temperatures corresponding to speeds exceeding approximately 3 km/s. For 
shock speeds below 3 km/s (T < 7000° K), the dissociation process is slower 
than the vibrational excitation rate, thereby allowing vibrational equilibrium 
to be achieved before dissociation becomes appreciable. At high shock speeds, 
Uqq > 6 km/s, additional processes like CO dissociation and ionization become 
more and more important. 

The flow-field profiles depend strongly on the rate of vibrational 
excitation (vibrational relaxation times). The effect of energy exchange 
between the modes is surprisingly small. The difference in density profiles 
for fast and slow energy exchange between the bending and the symmetric 
stretching modes is almost an order of magnitude less than the difference in 
profiles with vibrations frozen or in equilibrium. It appears that shock-wave 
profiles depend primarily on the rate at which the 'translational energy is fed 
into internal modes and not so strongly on the energy distribution among the 
modes. It follows that one cannot obtain the vibrational energy in different 
modes nor different relaxation times with any accuracy from experiments that 
measure overall quantities like density or translational temperature only; nor 
can such experiments determine the rate of energy transfer into different 
modes. This could be done only by a detailed spectroscopic study, observing 
vibrational transitions. If flow-field profiles for overall quantities such 
as the density or temperature are of primary interest, a model analogous to 
that used for diatomic gases should give satisfactory results, that is, only 
the sum of the vibrational energies over all modes is considered and only one 
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effective relaxation time is used. It must he noted, however, that the true 
picture is complicated hy the many rate processes involved in the kinetics of 
polyatomic molecules. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Aug. 17, 1 966 
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APPENDIX A 


SELECTION OF CHEMICAL REACTION RATES AND VIBRATIONAL RELAXATION TIMES 


Chemical Reaction Rates 

At the present time, knowledge of the chemical reaction rates behind 
shock waves in C0 2 is insufficient to describe the state of the gas adequately. 
If one assumes, as in this paper, that only C0 2 * CO, C, and 0 are present 
behind the shock wave there would be two chemical reactions, each having the 
possibility of four different collision partners. This would then entail a 
knowledge of eight different reaction rates. The existing experimental data 
are summarized below. 

C0 2 rates.- The dissociation rate of C0 2 has been studied by several 
investigators . Davies (ref. 11) studied the dissociation of C0 2 in an argon 
heat bath behind shock waves by monitoring the infrared and ultraviolet emis- 
sion. The infrared measurements appear to give the best results and indicate 
that the dissociation proceeds from higher vibrational states. The upper 
limit of the temperature range for these data is 11,000° K. 

Steinberg (ref. 12) also monitored the infrared emission behind shock 
waves and determined the dissociation rate. His measurements were at tempera- 
tures up to 6,000° K in argon and nitrogen heat baths and his rates compare 
well with those of Davies. He obtained an activation energy similar to that 
obtained by Davies, and again found that dissociation proceeds from excited 
states. Brabbs, Belles, and Zlatarich (ref. I3) measured the dissociation 
rate of C0 2 in an argon heat bath at temperatures around 2,700° K with results 
similar to those of Davies and Steinberg. For the analysis of this paper, the 
analytical expressions of Howe et al. (ref. l) were used. Howe T s expressions 
give dissociation rates below those of Davies and Steinberg by a factor of 5 
to 10 at 5,000° K and above the experimental data by about a factor of 2 at 
10,000° K. They were used since they allow a convenient comparison of the 
results of this paper with the results of Howe et al. 

CO rates . - The only published work on dissociation rates of CO applicable 
to this work is that of Davies (ref. Ik) who monitored both ultraviolet and 
infrared emission behind shock waves in CO in an argon bath. The temperature 
range of the experiments was from 9*500° to 12,000° K and the data were fitted 
using an activation energy considerably less than the dissociation energy. 
However, the magnitude of the dissociation rate of CO has little effect upon 
the analysis presented in this paper. Therefore, the analytical expression of 
Howe et al. (ref. l) was used. The following constants for the two reactions 
were used in the Arrhenius expression (eq. (26)) for their reaction rates: 


2k 



Reaction 


1 

2 


B f r ' E f r ' B c r ' E e r , 

cm 2 (X f r k cal g mo1 a c r k cal 

g mol sec °K g mol cm 3 °K a Cr g mol 

6 . 955X10 12 0.5 125.7497 0.50596X10 14 -2.9713 138.7454 

7 . 7206 X 10 12 .5 256.1742 . 34113 XIO 3 -.2366 257.4958 


Vibrational Relaxation Times 

Many investigators have measured the vibrational relaxation times for 
C0 2 at temperatures below 2000 ° K. Since there is considerable disagreement 
in the results, a critical evaluation is necessary to select the most reason- 
able data and to substantiate their extrapolation to the temperature range of 
this study. 

When only the experiments made in shock tubes are considered, all are 
characterized by measurements of some time -dependent gas-dynamic variable, 
usually density, and the data are reduced by inferring that the relaxation 
time must be that giving the observed profile. These measurements, in general, 
give only a gross relaxation time; the most important results are the tempera- 
ture dependence of the rates and the effect of impurities, mainly, water vapor. 
Smiley and Winkler (ref. 15) , Griffith, Brickl, and Blackman (ref. 1 6 ) and 
Daen and de Boer (ref. 17 ) found that the logarithm of the relaxation time, 
referred to one atmosphere pressure, was proportional to T “ 1/3 which agrees 
with theoretical predictions. There is about an order of magnitude difference 
between the above results, but Daen and de Boer, using carefully dried C0 2 , 
have shown that most of the disagreement can be resolved by allowing for the 
amount of water vapor in the test gas. 

Since C 0 2 has four normal modes, it is important to know if all of the 
modes have the same relaxation time. Griffith, Brickl, and Blackman (ref. 1 6 ) 
observed that the apparent downstream equilibrium density was lower than 
expected. This could be explained by the assumption that the stretching modes 
had not yet reached equilibrium, thus indicating the possibility of separate 
relaxation times for each mode. In a later note, Greenspan and Blackman 
(ref. 18) pointed out evidence for three separate relaxation times. Witteman, 
in an analytical study (ref. 19 ) j found two different relaxation times for the 
direct excitation of the bending mode and the symmetric stretching mode, the 
time of the latter mode being one order of magnitude greater than the former. 

Based on an evaluation of the above experiments and analytical studies, 
relaxation times at T = 1000° K were taken as t x = r 4 = 10r 2 and t 12 = t 2 
with the transition probability P 1_> ° the same for all modes. The data of 
references 15 and 1 6 were used to determine the transition probabilities and 
yielded the following expressions for vibrational relaxation times: 

r m p = K m e 37/Tl/3 (1 - e " 0m/T ) _1 (Al) 
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with 


Ki = K 4 = 2.2X10“ 


sec atm 


Ko = K 3 = 2.2X10' 


sec atm 


T , =t) - T' 1/a xlO' 4 [ei=.(Tp) - ei i (Tp ) ] 

Tl2P ' iSfri) 


sec atm 


(A2) 


The experimental data were obtained to 2000° K and contain a spread of about a 
factor of 2. The extrapolation of the data to higher temperatures is "based on 
the Landau -Teller theory. Since the validity of this theory at high tempera- 
tures is questionable, one cannot expect the results to he accurate to more 
than an order of magnitude . 
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APPENDIX B 


VIBRATION AND DISSOCIATION COUPLING MODELS 


The first model for vibration-dissociation coupling implies that the 
vibrational energy of each mode of a molecule is limited to D/4 independent 
of the total vibrational energy of the molecule. This requires that each v m 
be smaller than a number V m . Because of the finite quanta of vibrational 
energy, this postulate cannot be complied with exactly, but one must make an 
approximation in the form V m < D/(40m). This gives the numbers: Vi = 8, 

V 2 = V 3 = 18, V 4 =4. 

With the probability that the molecule has the vibrational energy e m 
(ref. 20) proportional to 


-Vi0i/Ti -v 2 0 2 /T2 -V303/T3 -V404/T4 

P = e e e e 


we obtain for the vibrational energy in mode m 
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and for the average vibrational energy lost by dissociation: 
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(Bl) 


(B2) 


(B3) 


(B4) 


The limiting value of the vibrational energy e for high temperatures follows 
from equation (B2) and is approximately: 
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lim 
T m ->°° 


-m 


(B5) 


m=i 


The exponential factor in the forward reaction rate, equation (28), becomes: 


' m 


s -E f /RT _ -D/RT 
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n 

m=i 




v m =° 
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(b6) 


m 




v m =o 


The product on the right side of equation (B 6) accounts for vibration- 
dissociation coupling. For frozen vibrations (T m -> 0), -► 0 which results 

in the effective activation energy Ef equal to the dissociation energy D. 
This result is expected, since dissociation is then proceeding from the ground 
state of vibration. The other limit of vibrational equilibrium (T m = T) gives 
T^ oo and a product term of equation (B 6) greater than 1. This means that 
the effective activation energy Ef must be smaller than the dissociation 
energy D. In the general case of vibrational relaxation (0 < T m < T), Ef 
is somewhere between these limiting values. The mathematical formulation of 
the second model for vibration-dissociation coupling takes into account the 
postulate that the total vibrational energy of a molecule is restricted to the 
dissociation energy. Since for the bending mode all functions can be written 
as functions of v 2 + v 3 , the summations can be reduced to triple sums. For 
simplicity, the new variable v s + v 3 is written as v 2 - With the upper limits 
of the summation as 


Vl = 1M1 


.Z a9 , 4 


0i 


(B7) 


Vo = 


(P/R) - V4.04 - V101 


02 


(B8) 


V 4 = 


D 

R0, 


(B9) 


one obtains for the vibrational energy in mode m 
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Vx=0 v 2 =o 



Vl=0 V 2 =0 


V4=0 



V4=0 



vgQa + V404 \ 

T 2 T4 / (v 2 + i)v m e m 


Vi9i , v g9g , Ja§a\ 

. T 1 T 2 ) (v 2 + 1) 


(Bio) 


The average vibrational energy is, respectively. 



Vi^O V2=0 V4—O 


(Bll) 


It must be pointed out that the vibrational energies now depend upon the 
temperature of all three different modes. The limiting value of § m for high 
temperatures can be calculated by Dirichlets* integral: 


lim 

r T ) ->oo 

- L m 
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The exponential factor in the forward reaction -rate coefficient (eq. (26)) 
for this model is 


RT 
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= e 



Vi v 2 ( Vi 9 r , V2Q2 . ^4^4 \ 
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For small changes of the vibrational energy the associated temperature 
changes can he calculated by solving a system of linear equations: 


d§ m = 


de m 


m' 


(BU) 


with Sem/dTjn from equation (BIO) . The vibrational energies e m are shown 
as functions of the vibrational temperatures T m in figure 6. 
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APPENDIX G 


EFFECT OF VIBRATIONAL EXCITATION ON THE FLOW-FIELD VARIABLES 
IMMEDIATELY BEHIND STRONG SHOCK WAVES IN C0 2 


In CO 2 the flow-field quantities Behind a shock wave depend much more on 
the amount of vibrational excitation than they do in air, since there are 
three different vibrational modes. 

The purpose of the following discussion is to show the values of the 
flow-field quantities of temperature, pressure, velocity, and density behind 
the shock front (x = 0 ) for the two limiting cases of frozen vibrations 
(frozen flow) and of vibrational excitation to the full classical value 
(equilibrium flow) . 

For high Mach numbers (M^ -> 00 ) y equation (12) yields: 


-Gl = 7 + 1 
P 7-1 

00 


(Cl) 


If f degrees of freedom are excited to the classical value at the 
beginning of the relaxation region, one obtains 
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f + 2 
f 


(C2) 


and 


f 


+ 1 


(C3) 


The degrees of freedom are counted as follows: translation, 3 ; rotation 

(linear molecules), 2; and vibration; 2 per mode. This gives 
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Frozen vibrations 


7 


= £ 

7 ' 



Vibrational equilibrium, 
diatomic molecules 
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_ 15. 
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Vibrational equilibrium, C0 2 


The strong shock approximations for the ratios of the flow-field quantities 
for vibrations frozen or in equilibrium for C0 2 with four vibrational modes 
(the values for a diatomic molecule are given in brackets for comparison) give: 
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P e f e + 1 

— - — = 2 ^ 

P f f f + I ^ ' 



P e f e (f-^ + l) 

p f = f f (f e + i) = 



f e (ff + i) 5 
ff(f e + X ) 2 


o.kd , 


[1-3 for diatomic gases] 

(ciO 

[0.75] 

(C5) 

[1.05] 

(06) 

[0.79] 

(ci) 


Behind the nonequilibrium relaxation zone, when the dissociation of C0 2 is 
complete, the flow-field quantities are the same for frozen and equilibrium 
vibrations, although the temperature and flow velocities at the beginning of 
the relaxation zone are more than twice as high for frozen flow than for 
equilibrium flow. 


It should be pointed out that for translational temperatures higher than 
about 16,000° K the assumption of vibrational equilibrium is no longer justi- 
fied, since the vibrational energy would be greater than the dissociation 
energy of C0 2 . 
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